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ABSTRACT : 

It is shown that by squeezing the vacuum fluctuations of the electromagnetic field the 
quantum fluctuations of the optical forces exerted on laser cooled two-level atoms, can be 
dramatically modified. Under certain conditions, this modification in concert with the enhanced 
average forces can lead to equilibrium temperatures below those attained under normal vacuum 
fluctuations. 

INTRODUCTION: 

Laser cooling of atoms in a quasi-resonant standing laser wave has been attracting 
considerable attention during the past few years^. Another exciting subject has been the 
modification of the statistical properties of the vacuum fluctuations of the electromagnetic field. 
Reduction of these fluctuations in one quadrature phase of the field by almost an order of 
magnitude has been already realized in the laboratory. It is well accepted that the minimum 
equilibrium temperature of laser cooled two level atoms is determined by the interaction with the 
vacuum fluctuations of the electromagnetic field.This raises the question whether the equilibrium 
temperature of two level atoms in squeezed vacuum can be lowered below the normal vacuum level 

and in particular below the so called “Doppler limit” of K^T= hT/2 for two level atoms. 

In the following, the physical origin of the optical forces in a standing laser wave is 
described and an intuitive model of the effects in a squeezed vacuum is offered, the modified force 
in squeezed vacuum is presented and compared to the force in a normal vacuum. In order to find 
the equilibrium temperature the modification of the fluctuations of these forces in squeezed vacuum 
is calculated. This calculation show, under certain conditions, a dramatic modification of these 
fluctuations relative to the normal vacuum state, it is found that, in an intense standing wave, the 
reduced fluctuations in concert with the enhanced average cooling force can lead to equilibrium 
temperatures below those obtained under normal vacuum fluctuations. Moreover, under certain 
ideal conditions even sub-Doppler temperatures may be reached. In the running wave case, 
however, the temperature can not be lowered below the normal vacuum level. In addition to being 
of potential use for laser cooling, these results offer an interesting glimpse into the quantum nature 
of the momentum exchanges between the atoms and the field. 
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A slowly moving atom (kv<r ) in a low intensity standing laser light wave experiences a 
velocity dependent force. This "radiation pressure" force is well understood in terms of absorption 
and spontaneous emission. As first envisioned by Hct'nsch and Schawlow^, the atom experiences 
an increased absorption of photons from the laser beam which is shifted closer to resonance due to 
the Doppler effect. This velocity dependent differential absorption can provide a cooling force for 
laser detuning to the red side of the atomic transition or a heating force for blue detuning. At high 
intensity, however, stimulated emission can change the sign of the force to a heating force at red 
detuning and to a cooling force at the blue side of resonance 3-4 . 



Figure 1: Probe absorption as a function of its detuning from a pump tuned 20T to the red of an atomic 
transition.A)at low pump intensity the probe sees higher absorption at positive detuning closer to the atomic 
transition B)In normal vacuum at high pump intensity the TWM process is induced leading to less absorption for 
frequency shifts closer to the atomic transition . C) At the same high intensity as in trace B but in squeezed vacuum 
the TWM process can change its lineshape leading to an additional cooling force.(Ref. 8 eq. 12 with N =0.1, 
M=0.33, 0= n and Q=8F)- 


This stimulated force(or "dipole force") has been explained within the framework of the 
dressed atom model 5 and equivalently as resulting from Two Wave Mixing 6 (TWM). The TWM 
resonance appears in pump-probe spectra as a dispersive lineshape (as a function of the probe’s 

detuning from the pump). This feature has a width of the excited state decay rate, T, and shows 
decreased absorption at probe detuning from the pump closer to the atomic transition (see figure 
lb). In this process the atom absorbs one photon from one wave and emits a photon into the 




counter- propagating wave, thus acquiring a momentum kick of 2h k. This process usually requires 
high laser intensity ; however, it has been shown to occur at lower intensity when the normal 
relation between the dipole decay rate I ~2 and the excited state decay T ( r2=0.5 T ) is modified by 

the inclusion of phase interrupting events ( T2>0.5 T ). This effect is due to the appearance of a 

TWM term at lower order in laser intensity proportional to (Tj/T- 0.5]P where P is the saturation 

parameter. This phenomenon is closely related to the dephasing induced extra resonances in Four 
Wave Mixing. These resonances , which originally have been studied by Bloembergen and co- 
workers, are induced whenever the normal decay rates of the the atom are modified. Their 
relevance to the stimulated force is discussed in more detail in reference 6. 

Armed with this insight into the connection between TWM and the stimulated force, it is 
instructive to find the effect of squeezing on the TWM process. Gardiner 7 has shown that in 
general squeezing the vacuum fluctuations results in two different decay rates for the two 
quadratures of the atomic dipole, one of which is larger and the other smaller than the normal T/2 
value in ordinary vacuum. Hence, after the decay of the fast decaying quadrature component of the 
atomic dipole one is left with the slowly decaying component which means that ^ can be much 
smaller than 0.5 T. One can therefore immediately see that the “extra resonant” TWM process can 

also be induced in squeezed vacuum. Calculation of the lineshape 8 ' 9 of this process shows that 
indeed the TWM lineshape becomes phase dependent and can even change to a "dispersive" 
lineshape with opposite sign (larger absorption closer to the atomic transition) as demonstrated in 
figure lc . This indicates that in squeezed vacuum the stimulated force can be induced at low laser 
intensity . Moreover, it can change sign to provide an additional cooling force instead of heating 
for red laser detuning from resonance. 


The physical system under investigation is a slowly moving two level atom (k v« F) in 

either a standing or a running wave (a motionless atom is considered in the fluctuations analysis). 
The atom is embedded in a broad band squeezed light, so that all of the modes coupled to the 

atoms are squeezed. The bandwidth of the squeezing is broad enough to appear to the atom as a 8 
correlated squeezed white noise. The correlation functions for the multi-mode squeezed field can 
then be written as 7 : 

<bt(t)b(t> = TN8( t - 1’ ) <b(t)b+(t’)> = T(N+l)5(t - 1’ ) Fn , 


<b(t)b(t’)>=<bt(t)bf(t’)>* = rMe('^ icoot + 2ik r ) 6( t - 1’ ) 
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Where b,bf are the operators defined in terms of the positive and the negative frequency parts of 
this field, N and M are the squeezing parameters, N is proportional to the number of photons in the 
squeezed vacuum, while M< N(N+1) signifies the amount of correlation between the sidebands 
and the equality maximum squeezing. In the following we will choose M to be real and positive . 

The Hamiltonian describing the interaction of the atom with the quantized multimode 
radiation field and a classical coherent field is given in the electric-dipole and rotating-wave 

approximation by Eq.2 

H =^ig > 0 a 22 + H rad - (pE^e' i(0t G + + op*E*e +i0X ) + h (o + b + ^o) 


where C0 q is the atomic resonance frequency, o 22 . 0=0, 2 and o+=o 2 , are the atomic operators , p 
is the atomic dipole moment , H 0 is the free Hamiltonian of the field, and E 0 is the amplitude of the 
coherent field. One can then find the master equation for an atom in squeezed vacuum and derive 

equations of motion for the atomic operators whose expectation value is given by: 


< o* 12 > = - y< Oj 2 > - rM< o 12 >* + Q< D > eq.3 

< D > = - T(2N +1)< D >+ T - 2[Q*< 0,2 >+Q<o, 2 >*] 

The average force can now be found by calculating, the expectation value of the atomic variables 
and subsequently the first order corrections due to the atomic motion 3 . This gives rise to the 
following expressions 11 for the expectation value of the optical forces acting on the atom in the 
running<F> and standing wave<F s > : 


(F> = 


fikTP 

2(2N+1+P) 


, , 2(2N+1)A , , x 
1h ( k v ) 

X(2N+1+P) 


Eq.4 


H =- 


ahAP 

0_(2N+1+P) 


4%Mcos( 2<»P + r 2 0_(2N+l-P) - 2xP 2 

1 (a v ) 

TxP_(2N+l+P) 


where P is the modified saturation parameter in squeezed vacuum given by P= 2|Q| 2 Q>J% , 
^e^pE/^i and the other quantities are defined by: r,=(2N+l)r , C>_ = 2N+1 - 2M cos(2<j>) , % 

= IyI 2 - r 2 M 2 , Y=r,/2 - iA , A = A + rMsin(24>), where A is the laser detuning from the atomic 
resonance. 
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It is instructive to examine the new expression of the force in the standing wave by 
comparing it to the force in ordinary vacuum. In this limit (N =M=0) the force is reduced to the 


well known expression of the force (ref.3 eq.18) given by: 


(f)= -ah A 


1+P 


1 


r 2 (l-P)- 2|y| 2 P2 

ra+P)# 


- a- v 


eq.5 


Note that in this limit the first term in the numerator of the velocity dependent part of eq.4 is zero 
while the other two terms are reduced to those of eq.5. The striking appearance of the additional 
term in squeezed vacuum is analogous to the result of ref.6. In this case, the introduction of 

classical phase noise results in the appearance of an extra term -4lyl2[T2/r- 0.5]P , (r2=r/2+r<j) 
where Ttj) is the rate of the phase interrupting events).This term can give the stimulated force at 
lower intensity when r 2 / r>0.5 as phase noise is added. 

Notice that in the case of T 2 / T< 0.5, this term can also be induced but with opposite 
sign. This is indeed the case with quadrature squeezing, which can result in either larger or smaller 
phase noise than the vacuum level. This in turn introduces two different decay rates for the two 
quadratures of the atomic dipole.One of these, T 2 x = T(N+M+0.5), is larger and the other, V 2 

X =T(N - M+0.5 ), is smaller than the normal T 2 =r/2 value. Therefore, the sign of the extra term in 

eq.2 4[lYl2-r 2 M 2 ]M cos(2<J>) P can be controlled by the relative phase <J> of the driving field with 

respect to the squeezed vacuum. Hence, the stimulated force can not only occur at lower laser 
intensity, it can change sign to provide an additional cooling force at red laser detuning from 
resonance.This modification of the force can be further correlated with the TWM lineshape which 
becomes strongly dependent on the laser phase (]) and can even change sign as indicated by our 
intuitive analysis . 

Other important modifications of the force in squeezed vacuum are described by the term, 
A+rMsin2<t>. This term gives rise to a force at zero detuning as well as strong variations of the 

force at small detuning ( A<TMsin2(])). These effects can be understood by noting that the 

dephasing induced lineshape of TWM at resonance is absorptive in normal vacuum, but it can be 
transformed to a dispersive lineshape in squeezed vacuum, giving rise to a force at resonance . In 

addition it has been shown that the TWM can have sub-natural linewidth at small detuning 710 . 
This indicates that one can obtain arbitrarily large cooling forces at small detuning as the number of 
photons in the squeezed vacuum N, and therefore the amount of squeezing , is increased. This 

can be understood by noting that V r < N - M+0.5 ) in the limit of N»1 becomes arbitrarily 
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small, r 2y = T / 8N. 

In the analytic solution shown above the force is calculated only to first order in velocity 
(i.e a linear velocity dependence is assumed), this is correct only for small velocities kv«r. 

Numerical solution of the OBE, however, can provide the full velocity dependence of the force. 
This solution is shown in figure 2 for ordinary (trace A)and squeezed vacuum (trace B) . This 
figure demonstrates that the stimulated force which gives a heating force in normal vacuum for 

velocities in the order of kv<T/2 (as expected from the TWM lineshape, fig. lb), can be 

transformed to a cooling force in squeezed vacuum. The dashed lines in the figure are the results of 
the spatially averaged analytic solution which show good agreement with the numerical solution at 
small velocities. 



Figure 2: The velocity dependence of the spatially averaged force, in normal vacuum trace A and squeezed vacuum 
trace B, obtained by numerical solution of the OBE. The dashed lines are the result of the analytic solution. The 

parameters used for this figure are: A=-3T, 0=1. 36 T, T=10 7 Hz and k=5890 A for both traces and N=l„ M=V2 
and 0 =0 for trace B. 

Figure 3 demonstrate the interesting dependence of the force on the driving laser phase ({) 

(using the analytic solution eq.2 with kv=T/2 ).This is shown for a constant number of photons in 

the squeezed vacuum, N=l, but for various values of the squeezing parameter M. Trace A plots the 
force for thermal light M=0 (i.e no correlation between the sidebands) with no variation on the 
phase, as expected. Traces B-D , however, show large variations of the force for increasing 

degree of squeezing up to the maximum value of M, (M 2 =N[N+1 ]). This dependence is due to the 
different amount of noise that the induced dipole sees at different quadrature phase. Figure 3 also 
shows that even a modest amount of squeezing induce large effects on the force. 





Figure 3: The spatially averaged force as a function of the laser phase <t> for increasing amount of correlation between 
the sidebands A)M=0 (thermal light no correlation), B) M=0.5 , C)M=1 and D)M=V2 (maximum squeezing). Other 
common parameters used are: A=-3T, Q=l. 5T and N=l. 

THE QUANTUM FLUCTUATIONS 

It was recognized by Einstein ^ as early as 1917 that the fluctuations of the optical force, 
originating from both spontaneous and induced absorption and emission processes, are important 
in determining the Maxwellian distribution of the atomic velocity in thermal equilibrium. A simple 
momentum diffusion process which comes readily to mind is due to the random direction of the 
spontaneous emission recoils . In addition to this geometrical source of fluctuations one should 
also consider the fluctuations in the number of photons emitted in a unit time . However, this 

process can have sub Poisson statistics, as shown by Mandel ^ in resonance fluorescence, and 
give rise to an anomalous contribution 14 w hich can decrease the momentum spread, as 
discussed by Cook 14. An additional momentum diffusion mechanism becomes dominant at high 
intensities in a standing wave due to fluctuations in the stimulated emission process between the 
counter-propagating waves 3,14,5 This photon exchange between the waves, results in a random 

transfer of 2/jk units of momentum to the atom. Finally as shown by Gordon and Ashkin an 
atom even in its ground state can have random recoils due to the zero point vacuum fluctuations. 

In the following it will be shown that the dynamics in squeezed vacuum modifies the fluctuations 
of all of these processes. 

We are now interested in finding the force fluctuations on a stationary atom which are 
given by the diffusion constant , 2Dp : 
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Eq.6 



dt [<F(0) F(t)> - <F> 2 ] 


Insertion of F(t) = — 10 VG + H.c. for the force (where G is the freely propagating field), using 


the correlation functions for squeezed vacuum, <G>=Q (since<b>=0) and the commutation relation 


aij oki=oii 8(j,k) for the atomic operators gives : 

(F(O)F(t)) = fi 2 [(at(0)o(t) + a(0)at(t)) l(VQ) 2 | - (at(0)at(t))(VQ) 2 - (a( 0)a(t))(Vft*) 2 


-K^k) 2 [ TN (ai i + a 2 2 ) + r(o 22 )] 8(t) 


Consider first the last terms which depend directly on the quantum fluctuations of the field. These 
terms in the limit of normal vacuum(N=0) can be transparently modeled ^ as the random 
instantaneous emissions of momenta h k at an average rate of r< 022 >- In the squeezed vacuum 


case, this effect is enhanced, by the absorption of squeezed photons rN<Oj j> and spontaneous 


emission rN<o 2 2 > due to the larger number of photons in the squeezed vacuum . 

In order to evaluate the remaining terms , which describe the effects of the interaction of the 
coherent field gradient with the atomic dipole fluctuations.we need to find the integral of the atomic 
dipole autocorrelation functions which after some algebra leads to the following expression for 
the diffusion in squeezed vacuum: 


2D^= Eb+fc 2 P : 


2^ 


2 (2N+1+P) 


[1 + H] 


fiVn— 


20- (2N+1+P) 


1 + Dj + D* + D*] 


Eq.8 


Where a=0 , p=k in a running wave; and a=-ktan(k x) , p=0 in a standing wave while the 
other terms are defined by : 0+ = 2N+1 + 2M cos (2(f)) 


Eb=fc 2 k 2 


f 


2N+P 


(2N+1-+?) 


+ 2N 


p 

1 

(2N+1+P) 2 

(2N+l)0- 


nr 

X . 


and the standing wave terms; 
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Di = 


4 P 

( 2 N+ 1 )X 

( 2 N+l+P) 2 d>- 0 + 

r 2 


A 2 if 

( 1 +- 1 ) 

ITi X 




8 P 2 


(2N+1+P) 2 0-<D + 


A" 


r 2 rr { <b- 


3 = 


4X 


p3 


rr 1 d>- 0 + ( 2 N+ 1 +P ) 2 


In the limit of normal vacuum ( N=M=0) Eq.8 corresponds to the the results of ref.3 eq. 30 


2Dp=ft 2 p 2 r 


2(l +P) 


1 +- 


(l+P> l If 


r 2 

1- — 


+ 


fi 2 cc 2 r- 


1 +- 


- 3 


2P 2 

H- — + 


2 ( 1 +p^ (i-fP) 2 lH 2 J (i+p ) 2 r 2 (i+p> 


J2 


P 3 


+ ft 2 k 2 r 


2(l4P) 


Let us first examine the terms in the running wave case by associating them with the normal 
vacuum limit. As we discussed previously spontaneous emission in squeezed vacuum, represented 
by the diffusion term Dq, gives rise to increased fluctuations as a consequence of the increased 
number of photons in this state. However, as can be clearly seen, the D ^ induced absorption 
contribution, even in the normal vacuum limit, can reduce the spread. This term has been shown 
to originate from sub Poisson statistics of the emitted photons. In fact in the normal vacuum limit 

Dj coincides exactly with Mandel’s Q parameter 13,14 : 

((An) 2 ) - (n) 

^ <") 

Q=0 indicates a variance of <n>U2 in the number of the emitted photons(i.e no correlation between 
the photons) and negative Q sub-Poisson statistics. Figure 4 plots Dj as a function of the laser 
detuning. In squeezed vacuum , with an appropriate phase between the laser and the squeezed 
vacuum, can reach a value of -1, whereas in normal vacuum the maximum effect gives D j = - 
3/4. This indicates that in squeezed vacuum the photons can be emitted in an orderly manner, thus 
eliminating this source of fluctuations. Unfortunately, in the traveling wave case, one can not take 
advantage of this phenomenon due to the increased spontaneous emission term Dq. This is 
demonstrated in figure 5a, which shows that the equilibrium temperature, given by 

KbT= Dp/(-3v<F>), increases with the amount of squeezing. 

We now turn our attention to the more complicated case of the standing wave, as in the 
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running wave case, we still have enhanced fluctuations due to the larger spontaneous emission 
term, D 0 . However, the much larger average cooling force, in conjunction with the smaller 
fluctuations, of the higher order terms, make it possible to reach temperatures lower than otherwise 
obtained in normal vacuum. Unfortunately, one can not take full advantage of both of these effects 



DETUNING ( r ) 


Figure 4: The deviation from Poisson statistics as a function of the laser detuning from resonance.The dashed line 
shows the maximum effect in normal vacuum , while the solid line indicates that almost no spread in the number of 
photons emitted in a unit time in squeezed vacuum can be achieved. The parameters used are Q -0.351" in the normal 

vacuum while Q =0.261", N=2, M^=N(N+1) and <J> = 0.5 tt in the squeezed vacuum. 



Figure 5: The equilibrium temperature as a function of the squeezing parameter N in: a) running coherent laser 
cooling wave with A=0.5T, i2=T , $ = 0.25 n. b) a standing wave with A=0, Q=5F , 4> = 0.6 n .The dashed line is 
the Doppler limit temperature in normal vacuum. 

at the same laser phase. Nevertheless, one can choose a particular phase which will minimize the 
temperature as the squeezing increases. Figure 5b demonstrates the lowering of the equilibrium 
temperature below the normal vacuum level (N=0) as the amount of squeezing increases. 


172 





Moreover, the temperature can be decreased to values slightly below the Doppler limit (which in 
normal vacuum is achieved at low intensity). Lowering the equilibrium temperature is not the only 
benefit of squeezing, in the above example, the average force becomes larger than the maximum 
value in the normal vacuum giving rise to a shorter equilibrium time. Notice, however, that high 
degree of squeezing is needed in order to reach sub-Doppler temperature, in addition , degradation 

from ideal squeezing (M^<N(N+1)) results in a temperatures higher than the Doppler limit. This 
unfortunately makes the experimental demonstration of this effect rather difficult. 

A few words are now in order to get some insight into the modification of the fluctuations 
in the standing wave case. We first discuss the various terms in the normal vacuum case. The most 
notable difference from the running wave is the appearance of higher order terms in P. These terms 

were interpreted as resulting from the fluctuations of the dipole force ^and become important 
at high laser intensity. In particular the P 2 term is the only term that does not saturate at high P. 
Hence, although one can use the stimulated force in normal vacuum at the blue side of resonance to 
give a very large average cooling force(with the advantage of very short equilibrium time) the large 

fluctuations make the equilibrium temperature much higher than the Doppler limit 

With regard to the modification of these processes in squeeze vacuum, we begin by 

comparing the Dj s term to its counterpart in normal vacuum. As discussed in the running wave 
case, this term can be interpreted as the deviation of the fluctuations from Poisson statistics. We 
found that the modified D x s in squeezed vacuum can reach values close to -1. However, the 
behavior in a standing wave is quite different from the running wave, as is the case in the normal 


vacuum . 

The next term, D 2 S , is not present in a running wave and we assume that it describes the 

fluctuations of the lower order stimulated force. In a previous publication it was shown that 
while the velocity dependent stimulated force does not usually occur at low intensity, the inclusion 
of phase interrupting events (which increase the dipole decay rate), gives rise to a stimulated force 
term at lower order in laser intensity. In order to show that this identification is correct and to get 
some insight into this term D 2 S was calculated in normal vacuum but with increased dipole decay 


17=172+17, where IT, is the rate of phase interrupting events, which gives: 


D^ = 


0 

2P' 2 


(1+P') 2 


2IT 


-1 


V A V 


LV 




+2 L 

r 


As can be seen an additional diffusion occurs as phase noise is added , T 2 >17 2. Note, that 

when r 2 <T/2, this term becomes negative. Analogously, D 2 S in squeezed vacuum may be 
associated with the fluctuations of the extra stimulated force in squeezed vacuum (the first term in 
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the average force in squeezed vacuum , <F S >. Moreover, the fact that one of the dipole’s 

quadrature components can decay at a rate smaller than 172 suggests that D 2 S in squeezed vacuum 
can become negative, as indeed is the case. 

Finally we turn our attention to the highest order term in P , D 3 S , which is associated with 
the fluctuations of the normal stimulated force.The modification of this diffusion term is critical for 
achieving lower temperature at high laser intensity where the stimulated force becomes dominant. 
Figure 6 shows the spatially averaged total diffusion constant 2Dp, in a high intensity standing 
wave, as a function of the laser phase, in normal and squeezed vacuum.This comparison 
demonstrates the dramatically reduced fluctuations in a high intensity standing wave under 
squeezed vacuum conditions. 



Figure 6: The diffusion constant in a high intensity standing wave as a function of the laser phase in the case of 
normal vacuum (dashed line) and squeezed vacuum(soIid line) .The parameters used are: A=0, Q=5T , N=5. 

As to the experimental verification of these interesting phenomena. Although 90% squeezing has 

already been achieved in the laboratory, it is important to note, that the calculation presented here is 
carried out with the assumption that the atom is embedded in squeezed vacuum. In practice, the 
output of present sources of squeezed light (degenerate parametric oscillators) can couple only to 

part of the 4k steradians enveloping the atom. A possible solution to this problem has been 
proposed by Gardiner who suggested coupling the squeezed modes to the atom in a micro cavity. 
The other important assumption here is that the spectrum of the squeezing is much broader than 

that of the atomic transition . Theories which include finite bandwidth of squeezing^ have shown 
that the essential features due to squeezing are preserved, even for a bandwidth of squeezing only a 
few times larger than the width of the atomic transition. It should also be noted that it will be 

interesting to develop the theory with bandwidth of squeezing larger than T but smaller than the 
Mollow sidebands separation . This can introduce the possibility of controlling all three of the 
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decay rates of the atom and therefore might reduce the problem the diffusion due to of the higher 
rate of spontaneous emission in a broad band squeezed light . 

In conclusion, this paper demonstrates a dramatic modification of the quantum fluctuations 
of the mechanical effects of light on atoms which are also embedded in squeezed vacuum. In the 
running wave case the temperature can not be lowered below the normal vacuum level. However, 
in conjunction with the modified average cooling force even sub-Doppler temperatures may be 
reached, under ideal conditions, for atoms cooled in a standing wave. These interesting results, in 
addition for being of potential use, offer some insight into the quantum statistics of the photon 
exchanges between the atom and the field under squeezed vacuum conditions. Further investigation 
of other schemes of cooling with squeezed light might also be beneficial 
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